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1  Introduction 


PDF  methods  have  proven  useful  in  modelling  turbulent  combustion,  primarily  be¬ 
cause  convection  and  complex  reactions  can  be  treated  without  the  need  for  turbulence¬ 
modelling  approximations.  Xu  &  Pope  [14,16]  and  Tang  et  al.  [17]  provide  recent  ex¬ 
amples  of  such  model  calculations  based  on  the  modelled  transport  equation  for  the 
joint  PDF  of  velocity,  turbulent  frequency  and  composition  (species  mass  fractions 
and  enthalpy). 

The  advantages  of  PDF  methods  in  dealing  with  chemical  reaction  and  convec¬ 
tion  are  preserved  irrespective  of  density  variation.  Since  the  density  variation  in  a 
typical  combustion  process  may  be  as  large  as  a  factor  of  seven,  including  variable- 
density  effects  in  PDF  methods  is  of  significance.  Conventionally,  the  strategy  of 
modelling  variable  density  flows  in  PDF  methods  is  similar  to  that  used  for  second- 
moment  closure  models  (SMCM):  models  are  developed  based  on  a  constant-density 
formulation,  and  then  the  same  models  are  used  for  variable-density  flows,  but  with 
density- weighted  statistics  in  place  of  Reynolds  averages.  However,  it  is  not  clear 
to  what  extent  this  strategy  is  legitimate.  First,  the  coefficients  in  these  models  are 
deduced  from  experimental  data  or  simple  turbulent  shear  flows  of  constant  density; 
second,  the  impact  of  density  variation  on  the  model  formulation  is  not  considered. 
However,  taking  account  of  density  variation  in  turbulence  models  is  not  an  easy 
task  due  to  the  lack  of  good  experimental  and  theoretical  understanding  of  density 
variation  effects. 

Second-moment  closure  models  for  variable-density  flows  have  been  advanced  by 
a  number  of  researches  [1,3,4,9,11].  Recently,  Lindstedt  and  Vaos  (1999)  [6]  im¬ 
proved  the  pressure-strain-rate  term,  specifically  the  isotropization  of  production 
model  (IPM),  by  considering  the  preferential  acceleration  effect. 

There  is  a  close  relationship  between  PDF  models  and  second  moment  closure 
models  which  has  previously  been  exploited  to  improve  the  modelling  in  both  ap¬ 
proaches  [2,8,13,18].  In  this  paper  we  develop  a  PDF  model,  incorporating  variable- 
density  effects,  which  is  consistent  with  the  SMCM  of  Lindstedt  &  Vaos  [6].  We  first 
consider  the  Reynolds  stress  model  and  deduce  a  consistent  PDF  model  for  velocity; 
we  then  deduce  a  consistent  model  for  the  turbulence  frequency. 


2  Second  Moment  Closure  Equations 

With  few  exceptions  (e.g.,  [9]),  moment  closure  models  solve  transport  equations  for 
density-weighted  quantities.  The  decomposition  as  used  in  [4]  is 

Ui  =  Ui  +  u ",  £  =  £  +  £",  (1) 

P=  (p)  +  p\  P={p)+p', 
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where  U  ,p,p  and  £  represent  the  fluid  velocity,  density,  pressure  and  a  scalar  vari¬ 
able  (e.g.,  mixture  fraction  in  nonpremixed  flames),  respectively;  and  ~  and  ()  de¬ 
note  density-weighted  (Favre)  means,  and  volume-weighted  (Reynolds)  means,  re¬ 
spectively. 


2.1  Exact  Reynolds  Stress  Equation 

To  simplify  the  notation,  we  use  Rij  to  denote  the  density-weighted  Reynolds  stress: 

Rij  =  u'fu'i  =  (pu"u'j)/(p).  (2) 

The  exact  Reynolds-stress  equation  for  variable  density  flow  can  be  written  [4,6]: 


=  (p)pH  +  %  +  <t>ij  ~  ( p)£ij , 


where 


J l  =  ±  +  fj± 
Dt  ~  at  ldxt' 


is  the  mean  convective  derivative; 

Tat  =  (fm1uX)/(p), 

is  the  triple  correlation; 

P  _  dU,  dU. i 


is  the  production; 
is  the  mean  pressure  gradient  term; 


=  -M>^r  -  «)a<P> 


dxj  N  j  dxi  ’ 


(3) 

(4) 

(5) 

(6) 

(7) 

(8) 


is  the  fluctuating  pressure  gradient  term;  and  ( p)£ij  is  the  viscous  dissipation  tensor. 
The  fluctuating  pressure  gradient  term  is  decomposed  as  [6] 


where 


and 


—  4>fj  +  <f>a  +  4>Ij, 

(9) 

$ ij  =  $ij  ~ 

(10) 

D  2  /  ,du'l\ 

dl)' 

(11) 

(12) 

It  may  be  observed  that  in  constant-density  flow  and  <pP  are  zero. 
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2.2  Modelled  Reynolds  Stress  Equation 

The  modelling  focuses  on  the  term  </>-.  Lindstedt  &;  Vans  discuss  the  terms  (j)f-  and 
4>Jj,  but  in  the  end  either  neglect  them  or  absorb  them  into  other  modelled  terms. 
Their  model  can  then  be  written 


faj  ~  {p)£ij  —  ( p){GuRtj  +  GjtRu  +  Coedij) 

+CAR(®ij  -  (13) 


where  the  tensor  Gij  and  the  constant  Co  =  2.1  arise  from  the  generalized  Langevin 
model  [2];  e  =  is  the  dissipation;  CAr  =  —3/10;  and 


dxi  ' 


(14) 


The  terms  in  Gij  and  Co  represent  the  constant-density  model,  whereas  the  term  in 
is  the  addition  to  account  for  variable-density  effects. 

For  future  reference  we  write  the  (partially)  modelled  Reynolds  stress  equation 
obtained  by  substituting  Eq.(13)  into  Eq.(3).  The  result  is 


/  ,\  DRjj  ,  dTjjj 
W  Dt  dxe 


(p)Pij  "b 

+  {p){GieRej  +  GjeRa  +  CoeSij ) 
+  CAr  . 


(15) 


2.3  Model  Equations  for  fc,  e  and  u 

The  exact  equation  for  the  Favre  mean  turbulent  kinetic  energy 

t  =  ^  =  i«X,  (16) 

is  obtained  by  taking  half  the  trace  of  Eq.(3).  Neglecting  the  term  \<j)u,  the  result  is: 

+  =(P)(P~e)  +  *-  (17) 

The  GLM  coefficient  Gy  satisfy  the  constraint 

GijRij  +  ^1  +  2^°)  £  ~  (13) 

so  that  half  the  trace  of  the  modelled  Reynolds-stress  equation  is  identical  to  Eq.(17). 
The  standard  model  equation  for  the  Favre  mean  dissipation  [4]  is 
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(19) 


<p)K  +  %  =  -  c^£)  + 

where  Te  is  a  transport  term,  and  the  values  of  the  constants  given  by  Lindstedt  h 
Vaos  [6]  are  C£\  =  1.44,  C&  —  1.92  and  Ce 3  =  1.2.  The  term  in  Ce3  is  of  course  zero 
in  constant-density  flow. 

We  also  consider  the  model  equation  for  the  Favre  mean  turbulence  frequency 
u  =  e/k  that  follows  from  the  above  equations: 

/  Jfa  =  iP[De_  _  (p)e  Dk 
KP}  Dt  k  Dt  k2  Dt 

„  /  p\ 

=  Tv+W&^Crt-Cuijj  +  Cua^j-,  (20) 

where  %,  arises  from  the  transport  terms  in  the  k  and  e  equations,  and  Cua  =  Csa  -  1, 
for  a  =  1,2,3,  i.e., 

Cwi  =  0.44,  CW2  =  0.92,  Cu3  =  0.2. 

3  PDF  Model  Equations 

As  usual  [7,19],  we  approach  PDF  modelling  from  the  Lagrangian  viewpoint  by  con¬ 
sidering  stochastic  models  for  the  position  X*(t),  velocity  U*(t),  turbulence  frequency 
uJ*{t)  and  composition  £*(£)  of  a  fluid  particle.  The  models  proposed  here  are: 

dX*(t)  =  U*(t)dt,  (21) 


+(Gij  +  Gfj)(U*  -  Uj)  dt  +  (CQe  +  E^dWi,  (22) 

and 

<2m*  (t)  =  -C3(u*  -u)ndt-  S„Du* (t)dt  +  [2C3C4u)Qw*]  1  dW.  (23) 

The  quantities  appearing  in  these  equations  are  defined  below  as  each  equation  is 
discussed  in  turn.  (A  specific  model  for  £*(£)  is  not  needed  in  the  subsequent  devel¬ 
opment.) 

3.1  Velocity  Model 

In  the  velocity  equation,  Gij  and  Co  are  the  same  GLM  coefficients  as  appear  in  the 
Reynolds  stress  model,  and  W(t)  is  an  isotropic  Wiener  process.  For  simplicity,  we 
do  not  consider  here  a  model  for  the  pressure  transport  [15]. 
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The  coefficient  Ca,  E  and  Gf-  are  to  be  determined:  if  they  are  set  to  zero,  then 
Eq.(22)  corresponds  to  the  usual  constant-density  GLM,  which  corresponds  to  Eq.  (15) 
at  the  Reynolds  stress  level  (for  constant-density  flow). 

The  coefficient  Ca  is  particularly  significant,  since  it  determines  how  the  mean 
pressure  gradient  accelerates  fluid  of  different  density,  (p*  denotes  the  density  of  the 
particle.)  For  Ca  =  0,  the  term  corresponds 


DU 

Dt 


-wv<p>- 


implying  that  all  fluid  elements  experience  the  same  acceleration,  independent  of  their 
density.  Conversely,  for  Ca  —  1,  the  term  corresponds  to 


DU 

Dt 


(25) 


implying  that  the  pressure  force  per  unit  volume  experienced  by  the  fluid  is  indepen¬ 
dent  of  its  density,  and  hence  less  dense  fluid  is  preferentially  accelerated.  At  first 
glance,  the  exact  Navier-Stokes  equations  suggest  that  Eq.(25)  is  the  appropriate 
model — but  this  ignores  the  phenomenon  of  “added  mass.”  As  the  less  dense  fluid 
accelerates,  so  also  does  the  more  dense  fluid  that  it  displaces.  Hence,  we  expect  the 
appropriate  value  of  Ca  to  be  between  0  and  1. 

From  the  stochastic  model  for  U*(t)  (Eq. (22)),  standard  techniques  [7j8]  can  be 
used  to  derive  the  corresponding  model  equations  for  the  mean  velocity  U  and  for 
the  Reynolds  stresses  Rij.  The  mean  momentum  equation  is 


dxi ' 


(26) 


This  equation  is  free  of  model  coefficients,  and  (apart  from  the  neglect  of  the  viscous 
term)  it  is  the  exact  Reynolds  equation  stemming  from  the  Navier-Stokes  equation. 
The  Reynolds  stress  equation  derived  from  Eq.(22)  is  : 


(p)  +  QgTijt  —  ( p)Pij  +  Ca$ij 

+{p)[(Ga  +  G*)Rej  +  (Gje  +  Gfe)Rei  +  (C0e  +  E)5ij].  (27) 


By  comparing  this  equation  to  the  model  of  Lindstedt  &  Vacs  (Eq.(15)),  we  see  that 
these  two  equations  are  identical  if  the  coefficients  Ca,  Gf-  and  E  are  chosen  to  satisfy 
the  constraint: 


+  Car  —  g  $<5^  —  Ca&ij  +  ( p)(E5ij  +  G£Rej  +  Gf£Rei).  (28) 

Clearly,  based  on  consistency  with  Lindstedt  &  Vaos’  model,  the  appropriate  choice 
of  Ca  is 
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(29) 


C„  =  1  +  CAr  = 

Then  the  consistency  condition  reduces  to 

§(1  -  C.ynSi,  =  (?>(£<%  +  G&R,j  +  Gf'Ru). 


(30) 


We  now  consider  three  ways  to  satisfy  Eq.(30)  exactly  or  approximately.  The 
simplest  specification  of  E  and  Gd  to  satisfy  Eq.(30)  is 

B=|(l-C.)*/0»),  Gts  =  0.  (31) 

This  is  not,  however,  an  acceptable  model  since  it  does  not  guarantee  that  the  stochas¬ 
tic  model  for  U *(t)  (Eq.(22))  is  realizable.  To  be  realizable,  the  coefficient  C0£  +  E 
must  be  non-negative,  which  is  not  assured  by  Eq.(31). 

The  second  specification  of  Gd  and  E  considered  is 

E  =  °'  G»  =  ^  >j>  V.  (32) 

where  R^1  denotes  the  inverse  of  the  Reynolds-stress  tensor  (such  that  R^Rjk  =  Sik)- 
As  may  be  directly  verified  by  substitution,  this  specification  exactly  satisfies  the 
constraint  equation,  Eq.(30). 

To  avoid  the  complexity  involved  in  the  inverse  R^1 ,  we  consider  the  third  speci¬ 
fication 

p  n  piA  _  Ca)$  p  (oo\ 

E  °’  °ij  2 (p)k  6ij' 

This  does  not  satisfy  the  constraint  equation,  Eq.(30),  but  it  does  satisfy  its  trace. 
Consequently,  compared  to  Lindstedt  &  Vans’  model,  Eq.(33)  leads  to  the  additional 
term  in  the  Reynolds-stress  equation 

SM  = . . .  (1  -  C„)»  (jle/fc  -  | .  (34) 


The  occurrence  of  k  in  the  denominator  of  Eq.(33),  and  similarly  R~l  in  Eq.(32), 
raises  questions  of  realizability  as  k  tends  to  zero.  However,  when  used  in  the  stochas¬ 
tic  model  for  U*(t),  both  specifications  are  realizable.  Using  Eq.(32),  for  example, 
the  relevant  coefficient  in  the  stochastic  model  (Eq.(22))  is 


g*(u;  -  v,)  = 


(i  -  cy* 


«5‘(t7  -  u,) 


(i  -  cy  dtp)  { (o 

3  (p)  dx , 


(35) 
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Each  term  in  braces  is  a  non-dimensional  quantity  which  is  bounded  as  k  tends  to 
zero  (or,  more  precisely,  as  the  smallest  Reynolds  normal  stress  tends  to  zero). 

To  summarize  the  principal  result:  the  stochastic  model  for  velocity  Eq.(22),  with 
Ca  =  1  +  Car  =  0.7,  E  =  0  and  Gy  given  by  Eq.(32)  is  realizable  and  corresponds  at 
the  Reynolds  stress  level  to  the  model  of  Lindstedt  &  Vaos. 

3.2  Turbulent  Frequency  Model 

The  stochastic  model  for  turbulent  frequency  (Eq.(23))  is  that  given  by  Van  Slooten, 
Jayesh  &  Pope  [12].  The  quantities  u*,u  and  Cl  are,  respectively,  the  particle  fre¬ 
quency,  the  (Favre  mean)  particle  frequency,  and  the  condition-mean  frequency  which 
is  designed  to  account  for  external  intermittency  and  whose  definition  is  given  in  [12]. 
For  fully  developed  turbulence,  Cl  equals  u,  and  in  general  the  dissipation  is  given 
by  e  =  kCl.  In  the  final  term  in  Eq.(23),  W(t)  is  a  Wiener  process,  independent  of 
that  in  the  velocity  model.  The  values  of  the  constants  give  in  [12]  are  C3  =  1.0  and 
C4  =  0.25. 

The  equation  for  the  mean  turbulent  frequency  deduced  from  the  stochastic  model 
(Eq.(23))  is  : 


<P> 


Du 


d 


\(p)u"u]  —  S^Clu. 


(36) 


Dt  dxi 

If  Cl  and  u  are  assumed  to  be  equal  (as  they  are  in  homogeneous  turbulence),  then  the 
source  term  in  this  equation  is  consistent  with  that  in  the  standard  model  (Eq.(20)), 
if  Sw  is  taken  to  be 

S„  =  C„1j-C„2-Cw3-^.  (37) 

The  first  two  terms  correspond  to  the  standard  constant-density  model  [12].  Hence, 
the  conclusion  is  that  the  stochastic  model  for  u*  is  consistent  with  Lindstedt  and 
Vaos’  model,  if  the  term  in  Cw3  is  included  in  the  specification  of  S^. 

The  occurence  of  k  in  the  denominator  of  Eq.(37)  again  raises  concerns  about 
realizability.  An  analysis  of  the  k  and  u  equations  shows  that  the  model  is  realizable, 
but  that  it  can  cause  k  and  e  to  become  zero  in  finite  time.  Because  this  is  difficult 
to  handle  in  numerical  implementations,  it  is  preferable  to  replace  $/(kCl)  in  Eq.(37) 
by  the  well-conditioned  approximation 

*  55  fcn  +  7|$|’  (38) 

where  7  is  a  small  positive  constant  (e.g.,  7  =  0.1). 


3.3  Scalar  Flux  Equation 

Given  a  stochastic  model  £*(£)  for  a  scalar  £,  a  modelled  equation  for  the  scalar  flux 
u"£"  can  be  deduced  from  the  models  for  £*  and  U*  [8]. 
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The  full  scalar  flux  equation  depends,  of  course,  on  the  mixing  model  contained 
in  the  model  £*(i).  But  we  note  here  just  the  terms  arising  from  variable-density 
modifications  to  the  velocity  equation.  With  Gf-  and  E  given  by  Eq.(32),  these 
contributions  are: 

<P> Wt “  ' ' '  “  +  |(1  -  CJflr1  (<)«]?"}  .  (39) 

4  Conclusions 

The  variable-density  turbulence  model  of  Lindstedt  &;  Vaos  [6]  has  been  extended  to 
the  level  of  a  PDF  model.  The  resulting  stochastic  model  for  velocity  (Eq.(22))  is  an 
extension  of  the  constant-density  Generalized  Langevin  Model.  The  coefficients  of  the 
additional  terms  in  this  model  equation  are  determined  by  requiring  consistency  with 
Lindstedt  &  Vaos’  model  (Eq.(29)  and  Eq.(32)).  The  stochastic  model  for  turbulence 
frequency  is  made  consistent  through  the  addition  of  a  source  term  (Eq.(37)). 

The  model  proposed  here  has  yet  to  be  implemented  and  evaluated  in  different 
circumstances  by  comparing  its  predictions  to  experimental  data;  and  it  is,  of  course, 
desirable  that  this  be  done.  But  nevertheless,  the  model’s  performance  can  be  ex¬ 
pected  to  be  reasonable,  in  view  of  its  close  connection  to  Lindstedt  &  Vaos’  model 
—  which  has  been  evaluated. 
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